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Abstract: We report new modeling and error reduction methods for 
differential-absorption optical-depth (DAOD) measurements of atmospheric 
constituents using direct-detection integrated-path differential-absorption 
lidars. Errors from laser frequency noise are quantified in terms of the line 
center fluctuation and spectral line shape of the laser pulses, revealing 
relationships verified experimentally. A significant DAOD bias is removed 
by introducing a correction factor. Errors from surface height and 
reflectance variations can be reduced to tolerable levels by incorporating 
altimetry knowledge and “log after averaging”, or by pointing the laser and 
receiver to a fixed surface spot during each wavelength cycle to shorten the 
time of “averaging before log”. 
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1. Introduction 

1.1 Overview 

Lidar remote sensing techniques are powerful tools for global measurement of atmospheric 
constituents and parameters [1, 2], The required precision and accuracy of such quantitative 
measurements have been increased to unprecedented levels for future airborne and space 
missions [3, 4]. To meet such stringent requirements, nadir-viewing, direct-detection, and 
pulsed integrated-path differential-absorption (IPDA) lidar techniques are being developed to 
measure the differential absorption optical depths (DAODs) of the target species in the 
column to the surface [5, 6]. From the DAOD and altimetry measurements and other ancillary 
data of the atmosphere (temperature, pressure, water vapor, etc.), the column number densities 
and dry mixing ratios of the target species can be retrieved [7, 8]. 

This paper addresses the following common issues with such DAOD measurements that 
have not been well addressed in previous publications: 1) error contribution from the optical 
frequency noise of a pulsed laser; 2) incorporation of a bias correction factor to reduce the 
averaging time for “log after averaging” (as named in [9]), and thereby minimize the impact 
of Etalon and atmospheric turbulence effects; 3) error contribution from the surface 
reflectance variations [10] when beam spots shift on the surface; 4) methods to reduce errors 
due to variations of surface height and reflectance; 5) estimation of the measurement precision 
of the transmitted laser pulse energy. 

We present a statistical model to address these issues, to provide insight and estimation of 
DAOD measurement errors and to guide lidar system design and development. The model 
addresses both types of DAOD measurement errors: the relative random error (RRE) and the 
relative systematic error (RSE) (i.e., relative bias). This paper is organized as follows. After 
introducing several representative IPDA lidar examples in the next subsection, we quantify 
the error contributions from the laser frequency noise in section 2. In section 3, additional 
detection noise sources and a bias correction factor are incorporated into the model. A 
numerical example of the model is presented in section 4. The impact of surface reflectance 
variation is evaluated in section 5. A new method to reduce errors due to surface height 
variation is described in section 6. More considerations for the lidar system are discussed in 
section 7. Some supporting details are provided in the Appendix. 

1.2 IPDA lidar examples 

The Active Sensing of C0 2 Emissions over Nights, Days, and Seasons (ASCENDS) mission 
[3] has been planned by NASA to measure the global distribution of carbon dioxide (COi) 
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mixing ratios (-390 ppm in average) to -1 ppm precision. A candidate IPDA lidar 
measurement approach for ASCENDS, being developed at NASA Goddard, allows 
simultaneous measurement of C0 2 and surface height in the same path [6, 7], As shown in 
Fig. 1, a pulsed laser is wavelength-stepped across a single C0 2 line at 1572.335 nm [11] to 
measure the DAOD. To make the measurements uniformly sensitive to concentrations in the 
lower troposphere, the online wavelengths (V 2 to v 7 ) are placed at the sides of the line where 
the absorption is mostly from C0 2 molecules in the lower troposphere [12]. The two-way 
transmittance of the C0 2 line from a 400-km orbit (Fig. 1, right) has a narrow linewidth (-3.5 
GHz) and steep slopes on the sides, making the measurements sensitive to laser frequency 
uncertainties. The offline points (v t and v 8 ) are placed in the adjacent window regions where 
the extremely low C0 2 absorption is essentially wavelength independent. To meet the mission 
end goal, the DAOD of the atmospheric C0 2 line needs to be measured with a RRE < 0.1% at 
multiple wavelengths across the absorption lines, and the variations of RSE need to be kept to 
a small fraction of the RRE. To reduce the RSE and RRE of the DAOD measurements, the 
detected photon count for each wavelength is averaged across multiple laser pulses before 
taking the logarithm of the on/off line ratio (i.e., “log after averaging”). This lidar approach 
will be used as a concrete example throughout this paper. 

The partial RRE contribution from the laser frequency noise needs to be < 0.03%, to allow 
more margins for other error sources. In this paper, we show that the standard deviation of the 
effective laser frequency noise at online points for C0 2 needs to be < 0.23 MHz to meet this 
target. To satisfy such stringent requirements, laser pulses at each fixed wavelength are 
generally carved by amplitude modulating a frequency stabilized continuous-wave (CW) laser 
[13] (Fig. 1, left). The 1-ps wide pulses need to be at least -100 ps apart, to completely clear 
the bottom -15 km of the atmosphere, which eliminates crosstalk from cloud scattering. 
Chirp-free amplitude modulation [14] is used so that the laser pulses retain the instantaneous 
frequency of the CW laser. 



Fig. 1. The laser transmitter provides the wavelength-stepped pulse train (left) to repeatedly 
measure at 8 points across the 1572.335 nm CO 2 absorption line (right). 


A similar C0 2 sounder space mission called A-SCOPE [4] has also been formulated by 
the European Space Agency. The A-SCOPE transmitter generates a pair of pulses, one at an 
online and the other at an offline wavelength, at an overall pulse pair repetition rate of 50 Hz, 
with a pulse energy of -50 mJ. Comprehensive lidar sensitivity and spectroscopic analyses for 
the A-SCOPE instrument operating at either 1.57 pm or 2.05 pm C0 2 bands have been 
recently published [5, 8, 15]. A similar approach has been adopted for a planned methane 
sounder space mission dubbed MERLIN [16]. 

2. Measurement errors from laser frequency noise 

Throughout this paper, x represents the ensemble average, <r( x) the standard deviation, and 
a 2 (x) the variance of x. 
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2.1 Treatment of laser frequency noise 

The instantaneous frequency v(t) of a single frequency CW laser can be modeled as an 
ergodic random process and the frequency noise is defined as S r (t) = v(t) — v. The electric 
field of a pulse carved from the CW laser can be represented by the analytic signal of the CW 
laser field multiplied by the real single pulse amplitude a(t-t Q ) centered at time t 0 , as given 
below 

u(t) = a(t-t o )exp(j[2xvt + 0(t)]}. (1) 

Here the instantaneous frequency is v{t) = v + (1 / 2jv)d(j)/dt. The single-sided energy spectral 
density (ESD) of u(t), denoted as L(y F ) as a function of the Fourier frequency v F , 
represents a short-term spectrum of the CW laser within a pulse duration At. The line-center 
frequency v c {t Q ) of L(v F ) is linked to the instantaneous frequency v(t) and ait) [17] as 

v e (f„) = L y F L N^y F ) dv F = v(t)h{t-t 0 )dt. (2) 

JO J —oo 

Here L N (v F ) = L(v F )J | L(y F )dv F and h{t) = a 2 {t)j | a 2 (t)dt are the normalized ESD of 

u(t) and normalized pulse intensity envelope, respectively. As to be seen in this subsection 
and in Appendix A, the laser frequency noise contributes to measurement errors through two 
complimentary factors, the line shape L(v F ) and the fluctuation of the line-center frequency 
v c (t 0 ). Both factors are jointly determined by the frequency noise and pulse amplitude a(t). 

As shown in Eq. (2), the fluctuation of the line-center frequency J rl (t 0 ) = v c (t 0 )-v is an 
average of the CW laser frequency noise 8 v (t) within a pulse duration At weighted by 
h(t-t 0 ). For a top-hat pulse shape, 8 vl (t 0 ) is simply a running average of 5 v {t) within At. 
The variance a 2 (S vl ) of S vl (t 0 ) can be obtained by integrating its power spectral density 
(PSD), which in turn can be derived from S Sv (f), the PSD of 8 v {t). Since 8 vl (t a ) is a 
correlation between 5 v (t) and h(t) as seen in Eq. (2), its PSD S Svl (f) is related to S Sv (f ) 
by 

S s M) = S Sv (f)\H(f)\\ (3) 

where H(f) is the Fourier transform of h(t). As the ESD of hit), \H(f)\ 2 serves as a 
window function to filter the frequency noise represented by S gv (f). For a top-hat 
h{t) = rect(t/Af)/Af, the corresponding window function becomes I //(/) l 2 = sinc 2 (/Ar). 

Fast frequency noise components above its effective bandwidth 1/A t are quickly averaged 
out as noise frequency / increases. As shown in subsection 2.5, this has been verified 
experimentally with the seed laser described in [13]. 

Formulation linking Liv F ) to the frequency noise S Sv (f) and the pulse amplitude a(t) is 
further derived and discussed in Appendix A. Interestingly, the remaining frequency noise 
that is filtered out by I H i f ) I 2 contributes to the bandwidth of the line shape Liv F ). 

2.2 Single pulse formulation 

The energy absorption of a single laser pulse through a thin layer of target molecules is given 
by 
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dL{y F ,l ) = -L(y F ,l)cr 0 (y F ,l)N(l)dl. 


(4) 


Here a 0 (y F ,l ) is the absorption cross-section of the molecules at Fourier frequency v F and 
path length /, and Nil) is the number density of the molecules at path length /. The laser pulse 
energy E(J) at a path length / can be obtained from Eq. (4) as 


£(Z) = £(0)exp[-r(v c> 0], 
r(v c ,Z) S j'a*(v e ,/W)d?'. 


(5) 


Here r(v c ,l) is the effective optical depth of the species of interest and a eff (v c , l) is the 
effective absorption cross-section defined as 


(• oo 

a eff (*T > 0 = J 0 (y F > l ) L l V (V'f . i)dv F 


T(v F ,l) = exp[-r 0 (v F ,/)], 
r 0 ( v F’ l ) = j' o cr 0 (v F ,l')N(l')dr. 


(• oo 

j 0 cr 0 (v F ,l)T(v F ,l)L(v F ,0)dv F 
j 0 T(v F ,l)L(v F ,0)dv F 


(6) 


Here t 0 (v f ,1) is the monochromatic optical depth of the species of interest, and other 
attenuation factors are assumed to be wavelength independent. 

To account for two-way absorption, the path length / is taken to be the accumulated 
distance that the laser pulse has traveled, running from 0 to 2 r G as the laser pulse travels from 
the spacecraft to the surface (outgoing path) and back to the spacecraft (return path). Here r G 
is the distance from the spacecraft to the surface. Our model uses the effective optical depth to 
include the line shape factor, which causes only a bias t(v c ,1)-t 0 (v c ,1 ) to the optical depth. 
As shown in Appendix A, this bias can be accurately corrected. When L(y F ) is much 
narrower than the spectral width of the target absorption line, this bias even becomes 
negligible, as shown in subsection 2.6 for the 1572.335 nm CO 2 line. When L{v F ) is broad, 
the one-way effective optical depth may no longer be the same for the outgoing and return 
paths so that r(v c ,2 r G ) = 2 r(v c ,r G ) may become invalid (see Appendix A for more details). 

Incorporating other attenuation factors [18], the detected laser pulse energy W s 
backscattered from the surface becomes 

w s =£A ex P[-r( v c’ 2r C )]’ (7) 

where E s is the transmitted laser pulse energy, A s = pT^ lm D / r G is the lump sum of 
attenuation factors excluding exp[-r(v c ,2 r G )], D is a instrument constant, p is the ground 

surface reflectance (in sr _1 ), and T atm is the one-way transmittance of the atmosphere 
excluding the target species. 

The lidar receiver can operate in either linear mode with an internal gain or photon 
counting mode. To cover both cases, the analog detector signal and noise are gain normalized 
and expressed in detected photon count K. For any type of radiation (including thermal light), 
the ensemble average and the variance of K and the optical energy W entering the detector 
area are related by [19, 20] 


K = aW, 

(8) 

a 2 (K) = F e K + a 2 a 2 (W). 

(9) 
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Here a = rj / (hv) is proportional to the detector quantum efficiency //, and h is the Planck 
constant. Note that only the shot noise contribution is multiplied by the detector’s excess 
noise factor F e . To estimate errors, we expand exp[— r(v c ,2r G )] in Eq. (7) around v to arrive 
at 


W s =E s A s exp[-T(v,2r G )][l-(dT/dv c )S vl (t 0 ) + bS vl \t 0 )], (10) 

where b = [(dr I dv c ) 2 -d 2 r I dv 2 ]l 2. From Eqs. (8) and (10), the ensemble mean of the 
detected signal photon count K s is 

K s = aW s =«£' s A s exp[-r(v,2r G )][l + /70- 2 (J„ 1 )], (11) 


where the variance a 2 (S vl ) of S vl (t 0 ) contributes a small fractional bias ba 2 (5 vl ) to K s . As 
illustrated in subsection 2.6, this bias is usually negligible in practice. Using Eqs. (9) and (10), 
the variance of K s is related to the laser line-center frequency noise as 

<r 2 (K s ) = F e Y s +cc 2 <J 2 (W s ) = F e Y s + (K s ) 2 ( dr / dv t ) 2 a 2 (S vl ) , (12) 


where the small contribution from the quadratic term bS v 2 (t n ) is neglected. 


2.3 Averaging across multiple pulses 

We now turn to estimate the two-way DAOD Ar = T(y m ,2r G )-T(Voff ,2r G ) between online 

and offline laser frequencies v on and v„ff . Up to section 6, r(y c ,2r G ) is assumed to be 
constant over the path and duration of the multiple pulses to be averaged. Although the 
measurement of E s is affected by the amplified spontaneous emission (ASE) in the laser 
amplifiers, the resulting measurement errors are negligible in practice and thus will not be 
included in the present model (see Appendix B for details). 

For each online or offline wavelength, the laser pulses are assumed to have the same h(t) 
and a period of t p . From Eq. (10), the following sum of normalized pulse energies (SNE) can 
be introduced to include the effect of the laser frequency noise 


ml ~hEU) exp 


[-r(v,2r G )l Y^A s (J) 
V «=o 


d t 

l-—Y(t) + bQ n (t) 
dv„ 


g[A0V t , 2 a + ^r p )], 

2 . 1.0 A (0-o 


(13) 


where the averaged laser line-center frequency noise S vn (t) across n pulses 
(i = 0, 1, 2, .., n - 1) is given by 

S n if) ^ § [ A s (i)S vl (t + ixt p )]. (14) 

L i= O A0« 

The mean and variance of S NEl are found to be 

=exp[-r(v,2r G )](^"; o 1 A i (o)(l + ^cr 2 (^ vl )), (15) 

^(S^) = (5^) 2 ( dr/dv c f a 2 ^). (16) 
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Superscripts (or subscripts) on and off will be used to indicate online and offline laser 
frequencies, respectively. The offline a 1 (S ^ :l ) can be neglected because of the negligible 
slope (dr/dv c ) off . 

Corresponding to S NE 1 . the sum of normalized photon count K s (SNK) is defined as 

S NK1 = E" = o K s (i)l E s (i) and its mean is found to be S NK1 = aS NEl . Here S NK 1 is named with 

the following scheme. S stands for sum, NK for normalized photon count, and the trailing 
digit 1 is used to distinguish S NK \ from S NK2 and S NK 3 to be introduced in sections 3 and 6, 

respectively. Using K s (i)K s (i') = a 2 W s (i)W s (i') [19], the variance of S NK \ is found to be 


c'OVi) = F'Snnki +a 2 a 2 (S NE1 ), 


J NKl' 

_^K s (i) 
mK1 hE s \i) 


— $nki I E sav — S K / E, 


2 

sav ‘ 


(17) 


Here S NNK i is named similarly, with NN indicating that K s is normalized twice. The integrated 
photon count S K and average transmitted pulse energy E sav are defined as S K 

and E sav = (^"_ 0 A(*))/ W > respectively. S NNKl is essentially S K normalized by E 2 av . The first 

term in a 2 (S NKl ) stems from the shot noise. The second term arises solely from the laser line- 

center frequency noise. 

From Eq. (15), At is related to S NKl by 


At = 


-In 


c°ff 
\ °NK 1 J 


+ ln(R A ) + ln (R b ), 


(18) 


where R A = A° n (i) / A'f (i) = 1 + S KA and 

R b = [1 + b on <j 2 (3 "" )] / [1 + b aff <j 2 (8"f )].For the rest of this paper, we assume the transmitted 
pulse energies are approximately the same for online and offline pulses so that 
S° K " = S'f exp(-Ar). At has been in effect estimated from the following “log after averaging” 
estimator [5-9] 


Ar„ 


-In 


?off 

v °nki y 


(19) 


To estimate errors, approximations ln(x) = lnfxj-(l/2)cr 2 (A:)/|xj (accurate to the third 
order) and cr 2 (ln(x)) = cr 2 (x)l(x^ (accurate to the second order) are used in our 


calculations. The bias b Al SNKl = At snk1 -At and variance of At Si 


are found to be 


= -\n(R A )-\n(R b ) + -^- 


C on nojf 

°NNKl °NNK\ 


C on 2 noff 

V °NK\ °NK\ 


+ 

dr ' 

2 

y 

V dv o) 


a 2 (.82), (20) 


2 /no, 

a ( Ar s«s-i)= 


a 2 (S 


off 

NK1 


7 off 


hA„. 

sf 


exp(Ar)]- 


dr 

dv„ 


v\82). ( 21 ) 
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As expected, the shot noise contribution to ct 2 (At snk1 ) is inversely proportional to the 

integrated photon count S K . However, it is important to use S NK \ rather than S K to estimate 
At, to remove the effect of variation in E s (i). The laser line-center frequency noise 
contribution to the DAOD RRE is [(dr/ dv c ) on / where cr(S"“ ) is the standard 

deviation of the averaged line-center frequency noise of the on-line laser (averaged across n 
pulses). It becomes clear that d°"(t) represents the effective frequency noise. Obviously, an 

on-line laser frequency bias Av„„ will result in a DAOD bias ~ (dr / dv c ) Av on [5]. This 
bias can be avoided by accurate calibration and thus is neglected throughout this paper. 

2.4 Frequency noise reduction from pulse averaging 



Frequency [Hz] 

Fig. 2. The nominal laser frequency noise PSD .S' , (/) = h f / f + h 0 (solid red) is 
suppressed to S g (f)\H(f)\ 2 (dotted blue) by averaging within a 1-ps top-hat pulse and 

further reduced to S Sv (f)\H{f)\~W an ff) (dashed black) by averaging across 1000 
identical pulses separated by a period of 1 ms. 


We now turn to link the PSD S Svn (f) of S vn (t) to S gv (f). From Eq. (14) we have 


s^(f) = s Sv (f)\H(ffw a M) 
W a Jf) = 


J «- 1 

Z A (0 ex P( ft p x 0 


( 22 ) 


Here the effect of multiple-pulse averaging on the PSD is an additional multiplicative window 
function W avn (f). The overall window function I H (f)\ 2 W avn (f ) is the ESD of the weighted 

pulse train function Z, "1[4 (f)h(t — ixt„ )] / Z"-o A O')- W OT (/) is a periodic function with a 
period equal to the pulse repetition rate. Using the orthogonal property of the phase shift 
factor e\p(j2/rft p xi), its averaged height is found to be 




(23) 


When the weighting factors A s (i) are the same, W avn (/) becomes sine 2 (nft )/ sine 2 (ft ) 
with an average height of 1 In. This is consistent with «-fold variance reduction when 
averaging n identical but uncorrelated random variables. As shown in section 5, < W avn > is 
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only slightly higher than 1 In when measured surface reflectance data are used for ,4/0, and 
will be taken to be I In for our numerical evaluations. 

The present model reveals relationships not available from previous models [5, 8]. Using 
Eq. (22), the effective frequency noise can be estimated in terms of S gv (f), as illustrated in 
Fig. 2 for a representative frequency stabilized seeder laser. The nominal S gv (f) can be 
approximated by S Sv (f) = h_ 1 /f + h 0 , where constant /; and h 0 quantify flicker noise and 
a nearly flat noise floor for /above 100 Hz, respectively. The periodical window function 
W a M) suppresses the frequency noise except in its narrow peaks at f = m w /t p 

(m w = 0, 1, 2, 3, ...). Frequency drifts with periods longer than the averaging time fall in the 
window peak at / = 0 and thus are not attenuated from this averaging. Since the PSD is 
nearly constant over the frequency range of all other window peaks (m w = 1, 2, 3, ...), the 
variance of the faster frequency noise components (excluding the slow drifts) is essentially 
reduced by a factor of n. This frequency noise reduction has also been verified 
experimentally, as shown in the next subsection. 

2.5 Experimental results 
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Fig. 3. (left) The standard deviation of the absolute frequency of the slave laser and frequency 
offset between the slave and master lasers measured within a gating time from 1 (js to 1 ms; 
(right) Both frequencies were also measured periodically within a 1-ps gating time with a 1-ms 
period and averaged across multiple measurements. The standard deviation of each averaged 
frequency, computed from 100 averaged samples, is plotted as a function of the number of 
measurements being averaged. 


Equations (3) and (22) have been verified experimentally with our seed laser described in 
[13]. Figure 3 (left) shows the standard deviation of the absolute frequency of the slave laser 
and frequency offset between the slave and master lasers measured within a gating time from 
1 (is to 1 ms. The absolute frequency noise was measured from the beatnote between two 
identical and independent laser units. The noise of the frequency offset was measured from 
the beatnote between the slave and master lasers. A fast frequency counter (HP 5371 A, 
Agilent Technologies, Inc.) was used to measure the beatnote frequencies. The measured 
noise decreases as the gating time increases, matching the reducing area under S gv (f)\H(f)\ 2 . 

Both frequencies were also measured periodically within a l-ps gating time with a 1-ms 
period and averaged across multiple measurements. The standard deviation of each averaged 
frequency, computed from 100 averaged samples, is plotted in Fig. 3 (right) as a function of 
n, the number of measurements being averaged. The variance of the averaged frequency 
decreases as n increases, matching the reducing area under S Sv (f)\H(f)\ 2 W m . n (f). 

2.6 Frequency stability requirements for ASCENDS 

Using Eqs. (5) and (6), the relative error (Ar-Ar 0 )/Ar 0 due to the finite line shape L(v F ) 
is numerically evaluated and shown in Fig. 4 (left) as functions of the effective bandwidth 1/A t 
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of a top-hat laser pulse shape for the C0 2 line. As detailed in Appendix A, the line shape 
L(v F ) is simply taken to be the ESD of the top-hat pulse shape. (Ar- Ar 0 ) / Ar 0 is smaller 
than a 0.02% target [5] and thus does not need to be corrected for all three online points (with 
frequency offsets + 0.5 GHz, + 1.1 GHz, and +1.7 GHz) when the pulsewidth At is longer 
than 0.5 ps. t is taken to be the same as r 0 for the C0 2 line in our numerical examples where 
At is assumed to be -1 ps. 



Fig. 4. Relative errors for atmospheric CO 2 DAOD measurements: (left) (Ar - Aro)/Aro due to 
the finite linewidth 1/A/ of a top-hat pulse shape (calculated from Eqs. (5) and (6)); (right) 
laser line-center frequency noise contributions to the RRE for G(S vn on ) = 0.1 MHz (dashed red, 
calculated from Eq. (21)), and RSE ln(R^)/Ar for a(<5 v i) = 10 MHz (dotted blue, from Eq. (18)) 
as functions of the online laser frequency offset from the absorption line center. The two-way 
optical depth r (solid green) is also plotted, with blue dots marking the same online and 
absorption line-center frequency points used on the left. 


Table 1. RRE Targets, Relevant Parameters, and Resulting Effective Laser Frequency 
Noise Requirement for Atmospheric CO 2 DAOD Measurement 


Requirement 

Value 

Parameter 

Value 

Budget for total RRE o(Ar) /Ar 

0.1% 

t(v c , 2 re) online (typical) 

~2.0 

Budget for partial RRE from laser 
frequency noise 

0.03% 

r(v c , 2 re) offline (15.6 GHz laser 
frequency offset) 

0.013 

Upper bound of G(S vn on ) 

(0.5 GHz laser frequency offset) 

0.23 MHz 

Sensitivity of RRE to <j(S vn ° n ) (0.5 
GHz frequency offset) 

0.13% /MHz 


The partial RRE [(dr/dv c ) on / Ar]cr(J°”) for <j( S °" ) = 0. 1 MHz, the partial RSE 
I ln(R 6 ) I /At, and t(v c , 2 r G ) are plotted in Fig. 4( right) for atmospheric C0 2 , as functions of 
the online laser frequency offset from the absorption line center. Table 1 summarizes the RRE 
targets, relevant parameters, and resulting effective laser frequency noise requirement for 
atmospheric C0 2 DAOD measurement. The frequency offset of the offline laser is set at 15.6 
GHz. To bound the RRE to < 0.03% for an online frequency > 0.5 GHz away from the 
absorption line center, the standard deviation cr(S°") of the effective frequency noise needs 
to be suppressed down to 0.23 MHz. This frequency stability requirement has been satisfied 
experimentally (see [13] and subsection 2.5). 

The partial RSE lln(R i )l/Ar for cr(S vl ) = 1 0 MHz is less than 0.002% at such online 
frequencies, and thus negligible. However, this bias and (At - Ar 0 ) / Ar 0 (Fig. 4, left) 
become significant near the absorption line center, where the excessive C0 2 absorption also 
results in higher RRE. It is thus undesirable to place an online frequency point near the 
absorption line center. Referring to Eq. (20), the last bias term (dT/ dv c ) 2 oll a 2 (8°”) / 2 is 
< 1(T 7 when [(dT/ dv c ) on I At]ct(S™) is bound to 0.03%, and thus negligible. The third bias 
term can be reduced to a negligible level with photon count integration. As shown in section 
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5, the remaining bias term lln(R A )l can be reduced to < 1.1x10 4 . Overall, the RSE 
b Ar snki / Ar can be reduced to < 0.02%. 

3. Incorporating additional detection noise sources 

We now incorporate additional detection noise sources arising from the background solar 
radiation, detector dark count, and receiver circuitry noise for the lidar receiver. The mean and 
variance of the photon count K„.„ produced by the background solar radiation is given by 

K n _ n = 2aN bgi B o At and cr 2 (K n _ n ) = F e K n n , where 2 N bgd is the PSD of the background 

solar radiation and B a the optical filter bandwidth. The equivalent detector dark count K d 
(normalized by the mean internal gain M e of the detector) has a mean and variance given by 

K d = AjAt and a 1 ( K d ) = F d K d , where 2 d is an equivalent dark count rate and F d an 
effective dark-count excess noise factor. The equivalent photon count K T arising from the 
receiver circuitry noise has a zero mean and a variance a 2 {K T ) given by 


c 2 (K t ) = I Si At / (M e e ) 2 , 

Is, = S si {f) sine 2 {f At) Atdf, 

J— oo 


(24) 


where S sj (f) is the PSD of the equivalent input noise current of the circuit and e the electron 
charge. Equation (24) is derived by the same steps leading to Eq. (40) for cr 2 (K sn ) (see 
Appendix B). The integration I Si can be approximated by S iU ((Y) when the detector electrical 
bandwidth B e is much larger than 1/Ar. 

The laser signal photon count K s is estimated as K\ = K tot - K ' bgd , where K tot is the total 

count K tot = K s + K nn + K d + K T measured within At, and K' bgd = K' n n + K' d + K' r is the 

background count within At scaled from a separate measurement taken between the pulse 
measurements. This background measurement can be taken in a longer duration f>At, where /> 
is a background integration time multiplier, to reduce the background variance. Obviously, 

K' s = K s . Referencing Eq. (9), the variance of K' s is found to be 

a 2 (K'J^FY + a 2 a 2 (Wj + \ d At, 

r (25) 

[ 2aN ^B 0 F e + F d A d + I Si / (M,e) J (1 + 1 1/3). 

The transmitted pulse energy E s can be estimated as E' s =E' tol - E' h?d where E tot is the 

total detected energy and the background energy E' bgd arises from the circuitry noise, dark 
current in the detector, and ASE in the laser amplifiers. The mean of E\ equals E s and the 
variance of E' s can be neglected (see Appendix B). The additional noise sources can be 

included by replacing S NKi with S NK1 = K' s (i) / E' s (i). Similarly, S K is replaced by 
S K , = 'y 1 " l K\ (i) . Obviously, S NK2 = S NKl . The normalized variance of S NK2 is found to be 


v 2 (S nk2 ) 



= F„ 


U NNK2 



f dr 


K dv cJ 


° 2 ($ vn ) 



(26) 


where S NNK2 and S NN2 are defined as S NNK2 = ^0 ’)/ E ' 2 (0 and S m2 = 

Here S NN2 is named similarly, with NN indicating double normalization. 
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As seen from Eq. (20), the “log after averaging” DAOD estimator in Eq. (19) leads to an 
additive bias term (1 / 2)F e (S°" NKI / S'f KI - S°f NK[ / S r f KI ) = (\l2)F e IS°f , and thus requires 
significant photon count integration ( S°“ ~10 4 ) to bring the bias term down to a tolerable level 

(~10“ 4 ). When using Snki , the additional background variance 2 hgd At in <j 2 (K' s ) results in an 

2 2 

additional DAOD bias term (1/ 2)A hgd At(S™ N2 / S°f K2 - S°f N2 / S'f K2 ). These two bias terms 
can be reduced significantly by incorporating the following bias correction factor 
C Ar snk 1 into the DAOD estimator given below 


^ T SNK2 — hi 


c . 


o°ff 
V °NK2 

( <?°ff 

°NNK2 
C off 2 
°NK2 


C A 


non 2 

^nk2 y 


Am 


( C"-'? 
°NN2 


n off 2 non 2 
\^NK2 ^NK2 y 


(27) 


As shown next, the sum of the two bias terms is reduced from ~(l/2 )F e I S°f, to 

— (1 / 2 )F e 2 / Sf! by the bias correction factor C Ar SNK2 . The residual bias - (1 / 2)F e 2 / Sf, 

becomes tolerable (~ 1 0 4 ) even when the integrated photon count .S'"" is as low as -100. Such 

a low photon count S°f can be obtained within a few pulses, or even within a single pulse 
without averaging. 

Using the bias-corrected DAOD estimator A t snk2 given by Eq. (27), the bias 
b M SNK2 = At snk2 - At and variance of A t SNK2 are found to be 


K.sm 2 = - hi(* A ) - ln(^ ) + (dT/dv e f on o 2 (SZ ) / 2 + b C AT _ SNK2 , 

1 1 — \ 2 3 — 3 ( 28 ) 

[l - exp(-2Ar)] - — F e nX bgd At\\ - exp(-3A r)] / S°f , 

O’ 2 (A t snk2 ) = o\ tDKT +(dT/dv c )l n o- 2 (#Z)’ 

2 (29) 

< DET = F. [1 ■ + exp(A r)] / Sf. + n\ gd At[l + exp(2A r)] / Sf, . 

It should be noted that all the parameters in the bias correction factor C Al SNK2 are 

measureable. For example, the background variance \ gd At can be measured when the laser 

pulses are off. However, it may not be feasible to measure some of the parameters (such as 
F e ) constantly during the flight. Consequently, long-term variations of such parameters may 
degrade the bias correction. 

To minimize the time- varying Etalon and atmospheric turbulence effects, it is desirable to 
reduce n, the number of pulses being averaged for At snk2 . To do this, the overall averaging 
time T m is divided into m slots (each is T tot / m long), and a A t snk2 (denoted as Ar SNK2 {k ) ) 
is evaluated within each slot ( k =0, 1, 2, .., m - I). At can then be estimated from the average 
of Ar SNK2 (k) given below 




,(*)■ 


(30) 
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The variance of A t snk2 av (excluding contribution from the slow laser frequency drifts) 
also decreases approximately by the number of participating pulses (i.e., nxm ), as in the case 
for A t snk2 . However, the averaging in Eq. (30) cannot effectively reduce the bias to At. 

4. Numerical estimation of measurement errors 



Two-Way C0 2 Optical Density 


— RRE (0.23MHz) 

— Shot noise 

— - Freqency noise (0.23MHz) 
Solar background 

— - Circuitry noise 

Dark count 

— - RRE target 
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Fig. 5. The RRE g(Atsnk 2 _civ ) /A r for atmospheric CO 2 DAOD measurement (solid black) as a 
function of the two-way optical depth, computed from Eqs. (29) and (30) using parameters 
listed in Table 2. The blue dots mark the three online frequency points used in Fig. 4. Also 
plotted are partial contributions to this RRE from the shot noise (solid grey), frequency noise 
(dashed red), solar background (dotted brown), receiver circuitry noise (dash-dotted green), 
and detector dark count (long-dashed blue). The standard deviation of the effective frequency 
noise o(<5 v „°”) averaged in 10 s is 0.23 MHz for the above calculations. When it is increased to 1 
MHz, the RRE (dotted black) rises above the 0.1% target (dashed black). 


Using realistic parameters listed in Table 2, At snk2 m / At, the RRE for atmospheric C0 2 
DAOD measurement, is computed and plotted as a function of x in Fig. 5. Partial 
contributions to the RRE from the signal shot noise ~ f’,[l +exp(Ar)]/ Sf , frequency noise, 
solar background, receiver circuitry noise, and detector dark count are also shown Fig. 5. 
When there are 3200 detected photons (in average) for each offline pulse, the RRE is < 0.1% 
for r between 0.66 and 3.2. This photon count can be achieved with E s = 4 mJ for an average 
p = 0.17. The detector specifications listed in Table 2 can be met with a state-of-the-art 
HgCdTe avalanche photodiode (APD) detector [21]. The detected solar radiation count K n . n is 
estimated for a zenith angle of 75°. Its noise contribution is larger than that from the receiver 
transimpedance amplifier (TIA). The noise from the dark count K d is negligible. It should be 
noted that the RRE increases sharply when the online point moves further away from the line 
center, confining the frequency offset of the online points to within +1.7 GHz from the C0 2 
absorption line center. The standard deviation of the effective frequency noise averaged in the 
10-s measurement time is taken to be 0.23 MHz for the above calculations. When it is 
increased to 1 MHz, the RRE rises above the 0.1% target (also shown in Fig. 5). 
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Table 2. Parameters Used to Evaluate the RRE of At as a Function of t 


Parameter 

Value 

Parameter 

Value 

Operating wavelength 

1572 nm 

Detector output dark current id 

1 nA 

Orbit altitude 

400 km 

Single-sided TIA equivalent input 
noise current 

1 pA/Hz 1 ' 2 

Telescope diameter 

1.5 m 

B e (electrical bandwidth) 

1.6 MHz 

Receiver field of view (full angle) 

150 prad 

B 0 (optical bandwidth) 

48 GHz 

Beam spot diameter on surface 

50 m 

Background solar photon count 
rate K n .JAt 

129 MHz 

T atm (one-way transmittance of 
atmosphere excluding CO 2 ) 

0.7 

P (background integration time 
multiplier) 

10 

p (surface reflectance) 

0.17 sr _1 

At (pulse duration) 

1 ps 

Receiver path transmittance 

0.51 

E s (transmitted pulse energy) 

4 mj 

A, (see Eq. (7)) 

1.48xl0~ 13 

K s (signal photon count) 

3200 

rj (detector quantum efficiency) 

70% 

T to t (overall averaging time) 

10 s 

M e (detector internal gain) 

400 

m (see Eq. (30)) 

10 

F e (excess noise factor) 

2 

effective n in T to t/m 

500 

Fd (dark-count excess noise factor) 

2 

averaged in T m 

0.23 MHz 


It is desirable to squeeze down RRE contributions from noise sources other than the signal 
shot noise. This will translate into reduction in required laser power and hence the cost. The 
laser peak power entering the lidar receiver needs to be high enough so that the signal shot 
noise contribution is well above the noise contributions from the receiver circuitry and the 
solar background count. This is the case for Goddard’s ASCENDS lidar where online points 
are placed at least 0.4 GHz away from the C0 2 line center with r < 2.4 . 

5. Impact of surface reflectance variation 


Table 3. Parameters Used for Evaluation of the Impact of Surface Reflectance Variation 


Parameter 

Value 

Parameter 

Value 

Speed of spacecraft 

7 km/s 

Pulse rate for each channel 

1 kHz 

Beam spot diameter on surface 

50 m 

Pulse duration 

1 ps 

Number of wavelength channels 

8 

Duration of averaging 

1 s 

Combined pulse rate 

8 kHz 

Path length within averaging time 

7 km 


As discussed earlier, the surface reflectance variation adds a bias lln(R A )lto the DAOD 
estimation. Since this bias varies from one path to another, it needs to be kept to a small 
fraction of the 0.1% RRE. Also, the PSD of the effective laser frequency noise 5 vn (t) is 
reduced by the periodical window function W avn whose average height < W [lm > depends on 
the variation of A s (i). Both factors are quantified here using measured surface reflectance data 
in [10] for realistic error estimation. 
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Fig. 6. (left) Relative surface reflectance: measured data (grey) and its running average over a 
length of 50 m (blue), courtesy of A. Amediek of Deutsches Zentrum fur Luft- und Raumfahrt 
(DLR); (right) ln(/? /4 ) calculated from surface reflectance data for the worst case. For each 

online wavelength channel, two ln(/? A ) are calculated from both offline sum 
for 1-s averaging time and the average is plotted (blue dot). An average of 10 consecutive 
ln(/? A ) values, each computed as above over 0.1 s averaging time, is also plotted (red 
square). 

The parameters used for the evaluation are listed in Table 3. The laser beam is assumed to 
travel at the same speed as the spacecraft. As shown in Fig. 1, the laser pulses are assumed to 
repeatedly cycle through wavelength channels 1, 2, 3, 4, 8, 7, 6, 5, 1, to minimize the 
separation between online (channels 2 to 6) and offline (channels 1 and 8) pulses. A set of 
surface reflectance data that has strong variations is used for this evaluation. The surface 
reflectance measurement was taken in southern Spain using ~10-m laser spot size and a step 
size of ~6 m [10]. To convert this data to the reflectance for our beam size, a 1-D running 
average is taken within our beam size (50 m) and the averaged reflectance is used as A s (i ) in 
our calculation. The raw and averaged reflectance data are plotted in Fig. 6 (left). For an 
averaging time of 1 s, there are n = 1000 pulses for each wavelength and the path length is 7 
km. In contrast, the separation between online and offline beam spots is only 2.1 m at most. 
Consequently, the ratio R A is very close to 1 . This point is confirmed by our calculation results 
shown in Fig. 6 (right). The 108-km total path length is divided into 15 sections (7 km each) 
and ln(R A ) is evaluated for each section. For each online wavelength channel, its starting 

pulse for sum y " !,A° n (i) comes in between two adjacent offline pulses, each serves as the 

starting pulse for its offline sum y ” _' n Af (i) . Each y ” 1 A°f (i) is used to calculate ln(R A ) 
and the average of the two is plotted in Fig. 6 (right) for the worst case where the bias 
I Inf W, ) I is < l.lxlCT 4 . The corresponding RSE is smaller than 0.02% if the two-way DAOD 
is kept > 0.55. 

I ln(/? A ) I is increased when the averaging time (hence n) is reduced. Figure 6 (right) also 
shows an average of 10 values of ln(/? A ), each computed as above over 0.1 s averaging time 
consecutively along the worst path. For some wavelength channels, this averaged bias 
I ln(/? A ) I becomes larger than 1CT 3 and thus unacceptable. To reduce the bias I ln(R A ) I, it is 
crucial to allow a fast pulse repletion rate (> 1 kHz for each wavelength) and a sufficient 
integration time (> 1 s), and to properly interleave online and offline pulses to minimize their 
separation. 

The average height < W avll > of W mn (f) is also evaluated for the 7-km averaging paths. 
The calculated < W am > is quite close to 1 In, varying between 1.01/n and 1.12/n. 
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6. Mitigating impact of surface height variation 

The DAOD errors due to surface height variation can be reduced to negligible levels if pulse- 
by-pulse knowledge of surface height is incorporated into the formulation. Let z 0 {i) be the 
surface height at laser beam spot i relative to a certain reference ground surface, and z(i) a 
measured value of z 0 (f). r G now represents the range from the spacecraft to the reference 
surface. The distance from the spacecraft to the actual surface is r G - z 0 (i). Since z 0 (i) could 
vary by tens of meters over the path of pulse averaging, r(r G - z 0 (;')) could fluctuate over a 
few tenths of a percent over the path, resulting in additional DAOD errors. To mitigate this 
problem, we estimate Ar(2 r G ) (i.e., the two-way DAOD from the spacecraft to the reference 

surface and back to the spacecraft) from S NK3 = K' s (i) /[A z (i)E' s (iJ], where 

A,(i) = exp[2a eJ/ {\',r G )N(r G )z(i)] is introduced to account for the fluctuation of the optical 
depth. The corresponding estimator At snk3 of At( 2 r G ) can be easily adapted from A t snk2 in 
Eq. (27) by simply replacing E'Ji) with A,{i)E' s (i) in the relevant terms in At snk , (i.e., 
S NK2 , S NNK2 , and S NN2 )■ As shown next, the variation of the surface height z 0 (i) does not affect 
such estimation of Ar(2 r G ). Nevertheless, the imperfect measurement of z 0 (i) still leads to 
additional DAOD errors. 

To simplify the treatment, we assume z(i) = z 0 (i) + Az + S,(i), where Az represents a 
slowly varying altimetry bias, and 5 r (i) a zero-mean random altimetry error. Az is treated as 
a constant during the n pulse averaging, and the SAi) values among different pulses are 
assumed to be uncorrelated and have the same variance cr 2 _ . Expanding 1 / A_ (/) in S NK2 as 

1 / A z (i) = exp(-2a eff NAz) / exp[2cr (# Nz 0 (i)] (l - 2a eff N8 z (i) + 2[a eff NS, O')] 2 ) , (3 1) 
the bias b Al SNK3 = A t snk3 - Ar(2 r G ) and variance of A t snk3 are found to be 

h —h 4 - h 

ty Ar_SNK 3 u At_SNK2 ^ u Atz * 

r - -|2 (32) 

Krz = 2(7 eff ( V °" > r G ) N 0b ) A ^ “ 2 O'™ , r G ) N ( r C ) (7 Z J C 1 " < W a m >) - 

& ( Ar SMT3) = Cr SNKl) + 17 Atz ' 

2 r - i 2 (33) 

O’ir. = [ 2(7 eff (^'on ,r Q )N (^ )(J z J < W am > . 

Elere we have neglected <7 cff {v off ,r G ) that is much smaller than cr eff (v 0 n,r G ). The additional 
DAOD bias /; Ar , arises from the altimetry bias Az as well as the variance a:. The additional 
DAOD variance al IZ scales with the altimetry variance al and < W am >, decreasing with 
~l/«. Taking atmospheric COt for example, an altimetry precision a =20m results in 
a i , ir (Von,r G )N(r G )( 20m) < 15xl0 4 , leading to a negligible standard deviation cr Arz < 1.3xl0 4 
(assuming n = 500), and a negligible DAOD bias ( < 1.4x 10 6 ) arising from cr 2 . On the other 
hand, the altimetry bias Az needs to be kept below 0.66 m to limit its DAOD bias 
contribution to < 10 4 . 
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7. Discussions 

7.1 Other error sources 

Other DAOD measurement error sources have also been studied previously [5-8, 15, 16]. An 
important one is the time-varying Etalon fringes in the lidar’s frequency response. It is crucial 
to minimize multi-path interference (MPI) along the optical train, and to shorten the time of 
“averaging before log”, to minimize the impact of the Etalon effect. Using multiple 
wavelength channels also helps to correct for the Etalon effect. The wavelength dependency 
of the lidar’s spectral response in the receiving path can be partially removed by passing a 
small fraction of the transmitted laser trough the receiving path and measuring the transmitted 
pulse energy E' s at the end of the path [15]. Another error in the measurement of K' s arises 
from the broadening of the returned signal pulses mainly due to the structured ground surface. 
The spectral crosstalk among multiple wavelength channels must be suppressed to maintain 
spectral purity for the measurement [5]. 

The contribution of atmospheric turbulence to the bias term I ln(J? A ) I is negligible because 
T atm essentially does not change during the time separation between an online pulse and either 
of its neighboring offline pulses (< 0.4 ms). The speckle noise can be also neglected in direct 
detection IPDA lidar instruments with a telescope diameter of ~1.5 m and laser beam spot 
size of -50 m on the surface [5, 7]. The Doppler effect can be either made negligible or 
corrected for practical nadir viewing IPDA lidars [15, 22]. The Doppler shift due to high 
speed cross-wind [15] can be made negligible by minimizing the off-nadir angle of the laser 
beam. The 2° off-nadir angle intended for A-SCOPE lidar to avoid detector saturation from 
the ice cloud glint [15] appears to be unnecessary because the reflectance of the ice (and 
snow) is particularly small for both 1.57 pm and 2.05 pm wavelength regions [23, 24]. The 
Doppler shift deduced from the spacecraft velocity components can also be compensated by 
dynamically shifting the frequency offset of the seed laser [13]. The laser beam divergence (~ 
+ 62.5 prad) in Table 2 results in a continuous spread of laser frequency of ± 278 kHz across 
the beam at 1570 nm [15]. Since this small spread is deterministic and symmetrical, it does 
not increase the variance of At and its bias contribution to At is negligible. 

7.2 Error reduction tradeoffs 

When the laser peak power entering the lidar receiver is high enough so that the signal shot 
noise contribution is well above the noise contributions from the receiver circuitry and the 
background count, the ratio rj / F e of the detector becomes a figure of merit for the DAOD 
measurement. For a fixed average power, further increasing the laser peak power by reducing 
the pulsewidth or the pulse repetition rate will no longer reduce o\ rDET effectively. On the 
other hand, it is desirable to use the maximum pulse repetition rate (-10 kHz for all 
wavelengths combined) to reduce the bias lln(/? A )l, and to use a wider pulsewidth (along 

with the high pulse repetition rate) to reduce <7 2 (S vn ) and the laser peak power. The laser 
peak power must be kept below the threshold of the undesirable stimulated Brillouin 
scattering (SBS) in the laser power amplifiers. The pulsewidth, repetition rate and laser peak 
power chosen for Goddard’s ASCENDS approach appear to be appropriate. Using a shorter 
pulse duration (< 100 ns) and a lower pulse repetition rate results in less frequency noise 
reduction, and may thus require a seed laser with a narrower frequency noise bandwidth and a 
lower noise floor. 

7.3 Shortening the time of averaging before log 

As shown earlier, all DAOD bias terms other than I \n(K A ) I can potentially be reduced to 
negligible levels by averaging across just a few pulses (or even without averaging) before log. 
It would be great if the same can be done for the bias term lln(J? A ) I, to reduce the time of 
averaging before log. This could be achieved by pointing the laser and receiver to a fixed 
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surface spot during each wavelength sweep cycle (< 1 ms). During such a short time, the slow 
spacecraft pointing jitter can be ignored. As long as the online/offline laser spots coincide 
within each wavelength sweep, slow pointing jitters do not increase lln(R A )l. The pointing 
can be achieved, for example, with piezoelectrically driven steering mirrors. While the 
spacecraft is traveling, the laser and receiver can be pointed to a fixed surface point for 
several online/offline cycles and then quickly redirected to the next fixed surface point. The 
mirror steering span can be reduced to only one or two beam spots to allow fast redirection. 
Even for a long pointing duration of 1 s, the fractional change in path length is less than 
4x10 5 and thus negligible. 

With this pointing scheme, the overall time of averaging before log can be significantly 
shortened to reduce measurement errors as mentioned earlier. On the other hand, using 
pointing mechanism may complicate the instrument design and pose a reliability challenge. 

8. Summary 

New modeling and error reduction methods are presented for DAOD measurements of 
atmospheric constituents using direct-detection IPDA lidars. Errors from the laser frequency 
noise are quantified in terms of the line center fluctuation and spectral line shape of the laser 
pulses, revealing new relationships that have been verified experimentally. Averaging across 
n pulses, the effective frequency noise is quantified by the multiplicative window function 
that filters the PSD of the instantaneous frequency noise. This in general leads to n-fold 
reduction in the variance of the effective frequency noise (excluding slow drifts with periods 
longer than the measurement averaging time). To bound the frequency noise contribution to 
the DAOD RRE to 0.03%, the standard deviation of the effective frequency noise needs to be 
< 0.23 MHz for ASCENDS’ C0 2 transmitter. 

The DAOD variance decreases approximately by the number of participating pulses. A 
significant DAOD bias can be essentially canceled out by incorporating a correction factor. 
The RRE and RSE of the DAOD due to surface height and reflectance variations can be 
reduced to tolerable levels by incorporating pulse-by-pulse altimetry knowledge and by “log 
after averaging”, or by pointing the laser and receiver to a fixed surface spot during each 
wavelength sweep to shorten the time of averaging before log. 

The predominant errors in the transmitted pulse energy measurement are found to be 
negligible in practice. Error reduction tradeoffs are investigated and illustrated with realistic 
calculations that take various detection noise sources into account. Other error sources are 
also summarized. 

Appendix 

A. Formulation linking L(v F ) to Sg v (f) and a(t) 

The bandwidth B L of L(v F ), as defined below, is also linked to v(f) and a(t ) [17] as 


Here A(f) is the Fourier transform of a{t) so that I A(f) I 2 is the ESD of a(t). 

As seen from Eq. (34), B 2 L is the sum of two components. The first component B LF 2 is the 
time domain “variance” of S v0 = v(t) -v c (f 0 ) weighted by h(t-t 0 ), and the second 
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component If 2 represents the bandwidth of I A( f ) I' arising from the variations of the 
amplitude envelope a{t). From Eq. (34), the ensemble average of B u 2 is found to be 


B LF 2 =[S v (t)f ~[S vl (t 0 )f =\y sv (f)[l-\H(f)\ 2 ]df. (35) 


This shows that B LF 2 arises from faster frequency noise components, excluding the slower 
components S Sv (f)\H(f)\ 2 that determine the line center fluctuation. Since the contributing 
fast frequency noise components are quite repeatable from one pulse duration to another, 
B lf 2 has little fluctuation around its ensemble mean B LF and thus causes little variance in 
b T (v c ,l) = T(v { ,l)-T t] (v c ,l). In practice, the frequency noise PSD S Sv (f ) cuts off eventually, 

resulting in a finite B LF 2 [17]. Since the slow frequency drifts do not contribute, B L 2 for a 
short pulse can be significantly smaller than that for the underlying CW laser. 

To facilitate further analysis, we rewrite L(v F ) as Lg{Av F ) = L{v c (t 0 ) + Av F ), as a 
function of the Fourier frequency offset Av F = v F — v c (t 0 ). Obviously, a broad L4(/)l 2 (due to, 
for example, a narrow pulsewidth) can still lead to a significant bias b T (v c ,l). Fortunately, the 
bias b T (v c ,l ) can be computed from Eqs. (5) and (6) using L 0 (A v F ) = L4(A v F )l 2 when the 
frequency noise contribution is negligible. Equations (5) and (6) thus allow accurate modeling 
even when L(v F ) is deterministically broadened to suppress SBS in the laser amplifiers. 
Since the amplitude modulation is deterministic, the bias b T {v c ,l) is deterministic and thus 
can be precisely calibrated and removed. Although L{v F ) for a top hat pulse shape has an 
infinite B^ , the relative error (Ar-Ar 0 )/Ar 0 can still be negligible as shown in Fig. 4 
{left). When L{v F ) is broad, the one-way effective optical depth may no longer be the same 
for the outgoing and return paths. This is due to the progressive distortion of the laser line 
shape L N (v F ,l) caused by the absorption of the species of interest. Since the laser line shape 
L N (v F ,l ) is no longer the same along the two paths, a c/j (v c ,l) may have two different values 

for each spatial position, one for the outgoing laser pulse and the other for the return laser 
pulse. 

The relationship between the CW laser line shape PSD and the frequency noise PSD 
S Sv (f ) has been established in previous studies [25-27]. The present work further links 
L{v f ), the laser pulse ESD, to S Sv {f) and the pulse shape a(t). L( v F ) is the Fourier 
transform of the following autocorrelation of u(t) 

/•co /• oo — 

^ u ( r o) = u*{t)u{t + t D )dt = a(t')a(t’ + t D )exp(j27rvt D )exp(jA<f>)dt' , (36) 

J —oo J —co 


where r' = r-f 0 and Ai/> = 0(t' + t o +t D )-t/>(f + t o ). The ensemble average of ^(Av^.) 
becomes 


L 0 (A V F ) = j t T a(t')a{t' + t D )exp ( j[A<j> - 2 xS vl ( t 0 )t D ])dt' 
Since A(j)-2nd vl {t Q )t D is a zero-mean Gaussian process, we have 


exp(-_/27rAv F r D )<* D . (37) 


exp ( j[ A (j) - 2 nS vl ( t 0 )t D ]) = exp 


-^(A0-27rS vl { t o )t D ) 2 


(38) 
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Using A<j . S = 


f *'+*<) + *D 
Jt'+t 0 


Sv(t)dt and taking similar steps as in [25, 26], we finally arrive at 


A> (A v F ) = f R 0 (t D ) exp(- JIkA v F t D )dt D , 

J — oo 

R 0 (t D ) = j ^a(t')a(t' + t D )exp — 2 tt 2 /‘ z , 2 J ^5 , ^(/)|sinc(/f 0 )-//(/)exp[y'2^/(r' + f D /2)]| 2 df dt'. 


(39) 


The phase factor exp[ j27if{t' + t D / 2)] in R 0 (t D ) stems from the time delay between A (f> 
(centered at t' + t D /2 + t 0 ) and 8 vl (t 0 ) (centered at f 0 ). Following the argument in [17] 
leading to Eq. (34), one arrives at B 2 L = -(1/ In) 1 (d 2 R 0 / dt 2 D ) lD=0 / /? 0 (0) that yields 

B 2 = B 2 ^ + B 2 f , the same result as given by Eqs. (34) and (35). The corresponding formula 
for CW laser [25] is a limiting case of Eq. (39) when a(t) is taken to be a constant. Similar to 
the CW formula, Eq. (39) in general leads to a non-analytical line shape L tl (A v F ) that can 
only be computed numerically. 

B. Impact of transmitted laser pulse energy measurement errors 

We now evaluate the DAOD errors arising from the transmitted laser pulse energy 
measurement errors. We assume that the laser amplifier ASE received by the detector has a 
uniform optical power spectral density of 2 N ASE ( N ASE for each polarization) within an 
optical bandwidth B o , and the laser is polarized and has a detected power of P s during At. For 
ease of analysis, we further assume that both the laser and ASE received by the detector are 
uniformly distributed among N MM spatial modes and the electrical-field product of any two 
different spatial modes does not contribute to the detector current (due to mode 
orthogonality). The detected photon rate then becomes A(t) = a[P s + P n _ n + P s _ n \ where 
P n ( t ) is the ASE power, P s n (t) corresponds to the product of the laser electrical field and 
the ASE electrical field (referred to as s-n beat term). aP s _ n (t ) has a zero mean and double- 
sided PSD S s _ n (f ) = 2a 2 1] N ASE / N mm within the detector electrical bandwidth B e [28]. The 
photon count variance arising from aP s _ n (t) is given by 


^\K s _ n ) = f 2 0tP s . n (t)dt] 2 = At J" Tri(-^-)r(g)d£ = I„_ n At , 

\J -At/2 ) J-cc At (40) 

I s . n ^r S s ^(f) sine 2 (fAt)Atdf. 

J —co 

Flere Y(c) = a 2 P t; n (t)l]_ n (t + g) is the autocorrelation of aP l; n it) and Trike) the triangular 
function. I s n can be approximated by S s _ n (0) = 2a 2 P s N ASE / N MM when B e »1/A t. On the 

other hand, P n _ n (t) produces a mean photon count K n n = a(2N ASI B o At) with a variance 
given by [29] 

<7\K n _ n ) = KZ(F e+ S c ), 

S c = KZj(2B o AtN MM ) = aN KS JN MM . 

To measure E s precisely, P s needs to be as high as -1 mW for a high signal to noise ratio. 
Such a high P s is also accompanied by a significant N ASI .. As a result, the variance cj 2 {E' s ) of 
E’ s arises predominately from the s-n beat term even though P n n (t) of the ASE produces a 
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photon count variance in great excess of the Poisson value. From Eq. (40), cr 2 (E' s ) is found to 
be cr 2 (E' s ) = 2 N ase E s / N mm . 

To account for the error 8E[ = E' s ~E s , we expand 1/ E’ s in S NK 3 (defined in section 6) as 
1 / E' s = (1/ £ s )[l- 8E[ / E s +(8E' s / E s ) 2 ]. Obviously, 8E' s values among different pulses are 
uncorrelated. The additional bias and variance of At S nk 3 arising from 8E[ are found to be 


N 
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(43) 


The additional variance <j A tE also scales with < W aln > 


1 In. In order to bound the 
corresponding RRE to < 1 0“ ' (assuming Ar ~ 2 , n = 1 and N MM = 1 ), the ratio N ASI . / E sav 
needs to be kept below 10 -8 . Consider a practical case where the seed laser oscillator is 
polarized with a peak power P seed and pulsewidth At = 1 ps at the input of a single mode 
preamplifier having a noise figure N F = 4. N ASE arises mainly from the ASE in the 
preamplifier that has an equivalent PSD of hvN F / 2 (in each polarization) at the input of the 
amplifier [28]. Using N ASE /E sav = hvN F / (2P seed At), the minimum seed laser peak power 
P seed at the input of the preamplifier is estimated to be -16 dBm, which can be easily 
satisfied. Similarly, the additional bias h Arl , can be neglected even when the online and 
offline pulse energies differ significantly. Consequently, both the bias b ArE and variance 
°ArE > arising from the laser- ASE beat noise, are neglected in the present model. 


C. Noise from background solar radiation 

At the lidar receiver, the ASE from the laser transmitter becomes negligible even compared to 
the weak background solar radiation that has a low optical power spectral density 2 N b d . 

Equations (40) and (41) can still be applied to the background sun light, provided that N MM is 
interpreted as the number of the received speckles and N ASE is replaced by N bd . Since 
aN hgd / N mm <k 1, K n n from the background sun light becomes a Poisson process (for F e = 
1). The s-n beat term variance cr 2 (K sn ) = 2aN bgd K s /N MM is much smaller than the signal 
shot noise F e K s and thus can be neglected. For the background sun light, the shot noise 
variance cr 2 (K n _ n ) = F e K n _ n could be reduced to a 2 (K s n ) if the optical filter bandwidth B a 
could be reduced to aP s / [(l + l/ P)N MM F ]. In practice, the optical filter bandwidth is much 
wider so that ct 2 (K n n ) is much larger than c 2 (K s n ) . 
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